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This paper introduces an experimental probe of the sterile neutrino with a novel, high-intensity 
source of electron antineutrinos from the production and subsequent decay of 8 Li. When paired with 
an existing ~1 kton scintillator-based detector, this (_E„}=6.4 MeV source opens a wide range of 
possible searches for beyond standard model physics via studies of the inverse beta decay interaction 
v e +p — > e + + 7i. In particular, the experimental design described here has unprecedented sensitivity 
to v e disappearance at Am 2 ~ 1 eV 2 and features the ability to distinguish between the existence 
of zero, one, and two sterile neutrinos. 

PACS numbers: 14.60.Pq, 14.60.St 



The beta decay-at-rest of 8 Li produces an isotropic 
electron antineutrino flux with an average energy of 
6.4 MeV. An underground liquid scintillator based de- 
tector can be used to detect these antineutrinos via the 
inverse beta decay (IBD) process v e + p — > e + + n. The 
antineutrino rate and energy, peaking at 9 MeV, can be 
fully reconstructed by the detector. Precise energy and 
vertex reconstruction opens the possibility of searching 
for antineutrino disappearance due to oscillations, which, 
in the simplest two-neutrino form, has the probability 

P = 1 - sin 2 26> sin 2 [1.27Am 2 (L/£)] , (1) 

where 8 is the disappearance mixing angle; Am 2 (eV 2 ) 
is the squared mass splitting; L is the distance (in me- 
ters) from the antineutrino source to the detector; and 
E (MeV) is the antineutrino energy. This probability 
is maximized in the range of Am 2 ~ E/L. An exist- 
ing large scintillator-based antineutrino detector with a 
diameter of O(10 m), when combined with an 8 Li iso- 
tope decay-at-rest source, is sensitive to oscillations at 
Am 2 ~ 1 eV 2 . This is an oscillation region of high in- 
terest due to anomalies that have been observed in the 
data from LSND [1], MiniBooNE [2], short-baseline re- 
actor studies [3J, and gallium source calibration runs g]. 
These anomalies are often interpreted as being due to 
sterile neutrinos [SHE] and have motivated the develop- 
ment of the IsoDAR (Isotope Decay- At-Rest) concept. 

IsoDAR-style sources have been considered before [9]- 
lllj . The design presented here, consisting of an ion 
source, cyclotron, and target, is the first with a suffi- 
ciently high antineutrino flux to address the existence 
of one or more sterile neutrinos. The 9 Be target, used 



mainly as an intense source of beam-induced neutrons, is 
surrounded by a 7 Li sleeve. When the target and sleeve 
combination is placed next to a kiloton-scale scintillator 
detector (e.g., KamLAND [T5J, SNO+ [13], or Borexino 
[14]). the large antineutrino flux from 8 Li beta decay can 
result in the collection of over 8 x 10 5 IBD interactions in 
a five year run. Such events allow a definitive search for 
antineutrino oscillations with the added ability to distin- 
guish between models with one and two sterile neutrinos. 
A sample of more than 7200 ?y e -electron scatters is also 
accumulated during this time and can be used as a sen- 
sitive electroweak probe. 

The charged particle beam, used for electron antineu- 
trino production, originates with a 60 MeV/amu cy- 
clotron accelerating 5 mA of ions. The design of this 
compact cyclotron [15] is ongoing and is envisaged as the 
injector for the accelerator system of the DAE<5ALUS 
physics program [T6l [17] . The IsoDAR design calls for 
about a factor of six increase in intensity compared to 
compact cyclotrons used in the medical isotope industry. 
Current and future medical isotope machines accelerate 
protons to 60 MeV and beyond with average intensities 
of 0.8- 1.6 mA PUE*]. 

In our design, a 5 mA H^" beam is injected at 70 keV 
(35 keV/amu) via a spiral inflector. Existing ion sources 
are sufficient to supply the beam required [20] . In the 
low-energy regime, space charge is crucial in modeling 
the beam dynamics correctly. The generalized perveance 
of a non-neutral beam, K = ql /(2ireom(3 3r Y 3 ), parame- 
terizes the strength of the space charge effect [21]. We 
conclude that the K for this machine is of the same order 
as existing proton machine designs based on ss2 mA and 
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an injection energy of 30 keV [221 [23]. The proposed high 
power is therefore feasible with regard to space charge 
issues. 

The accelerator described is a continuous-wave source 
with a 90% duty cycle to allow for machine maintenance. 
In consideration of target cooling and degradation with 
600 kW of beam power, we require a uniform beam dis- 
tributed across most of the 20 cm diameter target with 
a sharp cutoff at the edges. Third-order focussing ele- 
ments in the extraction beam line are able to convert 
the Gaussian-like beam distribution into a nearly uni- 
form one [23] and hence create the necessary condition 
on the target. 

The 60 MeV proton beam impinges on a cylindrical 
9 Be target that is 20 cm in diameter and 20 cm long. 
The primary purpose of this target is to provide a copi- 
ous source of neutrons. Neutrons exiting the target are 
moderated and multiplied by a surrounding 5 cm thick 
region of D2O, which also provides target cooling. Sec- 
ondary neutrons enter a 150 cm long, 200 cm outer di- 
ameter cylindrical sleeve of solid lithium enveloping the 
target and D2O layer. The target is embedded 40 cm 
into the upstream face of this volume; a window allows 
the beam to reach the target. The sleeve is composed 
of isotopically enriched lithium, 99.99% 7 Li compared to 
the natural abundance of 92.4%. The isotopically pure 
material is widely used in the nuclear industry and is 
available from a number of sources. The isotope 8 Li is 
formed by thermal neutron capture on 7 Li and to a lesser 
extent by primary proton interactions in the 9 Be target. 
For enhanced production, the sleeve is surrounded by a 
volume of graphite and steel acting as a neutron reflector 
and shield. The volume extends 2.9 m in the direction 
of the detector. Isotope creation in the shielding is neg- 
ligible. Figure [T] displays the target and sleeve geometry, 
and Table [T] summarizes the experimental parameters. 
We note that the geometry of the design is similar to 
that described in Ref. [TU] , 
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FIG. 1: A schematic of the IsoDAR target and surrounding 
volumes. The dots represent 8 Li (p e ) creation points, ob- 
tained with 10 5 60 MeV protons on target simulated. The 
neutron reflector, shielding, and detector are not shown. 



Accelerator 


60 MeV/amu of H+ 


Current 


10 mA of protons on target 


Power 


600 kW 


Duty cycle 


90% 


Run period 


5 years (4.5 years live time) 


Target 


9 Bc surrounded by Y Li (99.99%) 


V source 


s Li /3 decay ({E„)=6A MeV) 


Fe/lOOO protons 


14.6 


V e flux 


1.29X10 23 V e 


Detector 


KamLAND 


Fiducial mass 


897 tons 


Target face to detector center 


16 m 


Detection efficiency 


92% 


Vertex resolution 


12 cm/y/E (MeV) 


Energy resolution 


6A%/y/E (MeV) 


Prompt energy threshold 


3 MeV 


1BD event total 


8.2xl0 b 


I7 e -electron event total 


7200 



TABLE I: The relevant experimental parameters used in this 
study. 
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FIG. 2: The expected antineutrino flux and detected event 
rate in the experimental configuration considered. The an- 
tineutrino flux mean energy from 8 Li is 6.4 MeV. There are 
8.2xl0 5 reconstructed events expected from the 1.29xl0 23 v e 
created in the target and sleeve in five years. 



We determine isotope production rates using a 
GEANT4 simulation [35]. Due to its vast range of ap- 
plications, GEANT4 provides an extensive set of data- 
based, parametrized, and theory-driven hadronic mod- 
els, each one specializing in different types of interactions 
within a specified range of energy. The QGSP-BIC-HP 
physics package was chosen for this particular applica- 
tion. The applicable physics model is the pre-compound 
nuclear one which is invoked by the Binary Cascade sim- 
ulation. Simulated hadronic processes include elastic 
scattering, inelastic scattering, neutron capture, neutron 
fission, lepton-nuclear interactions, capture-at-rest, and 
charge exchange. For neutron energies below 20 MeV, 
the high-precision package uses the ENDF/B-VII data 
library [2"6] . 

Although all isotopes are considered in this analysis, 
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the induced 8 Li source in the sleeve dominates the an- 
tineutrino flux. The simulation yields 14.6 8 Li isotopes 
for every 1000 protons (60 MeV) on target. Approxi- 
mately 10% of all 8 Li is produced inside the target; the 
rest is produced in the sleeve. Neutrinos and antineu- 
trinos from other unstable isotopes are produced at a 
comparatively negligible rate. Over a five year run pe- 
riod and with a 90% duty cycle, 1.29xl0 23 antineutri- 
nos from 8 Li are created. IsoDAR's nominal oscillation 
analysis is done in terms of "shape-only" in L/E and is 
therefore independent of the absolute flux normalization. 
However, a "rate+shape" analysis using an absolute flux 
normalization uncertainty of 5% is also considered in this 
study. 

The IsoDAR antincutrino source is paired with an ex- 
isting underground scintillator-based detector for char- 
acterizing the antineutrino flux as a function of distance 
and energy. As can be seen in Eq. [TJ a baseline of 
L ~ 10 m is appropriate as a probe of the Am 2 ~ 1 eV 2 
region given the antineutrino spectrum shown in Fig. [2] 
We assume the face of the target sits 16 m from the 
center of a KamLAND-inspired 897 ton detector when 
calculating rates and oscillation sensitivity. 

The v e events are detected through the IBD interac- 
tion, which is unique in several ways. IBD has a com- 
paratively high cross section (~ 3 x 10~ 42 cm 2 ) which 
is known to < 1% across the relevant energy range [27] , 
The distinct delayed coincidence signal of prompt light 
from the positron, followed by a 2.2 MeV gamma from the 
neutron capture on a proton, enables a low background 
rate expectation. The antineutrino energy can be fully 
reconstructed on an event-by-event basis using the visi- 
ble energy of the prompt positron signal in the detector: 
Ev e = E e + + 0.78 MeV. A prompt energy threshold of 
3 MeV is employed here. In a five year run, 8.2 x 10 5 
reconstructed signal events are expected with 92% de- 
tection efficiency [2 8) and in the absence of antineutrino 
oscillations. 

The unique delayed coincidence signal makes reac- 
tor antineutrinos the only significant background in this 
analysis. The reactor antineutrino event rate at Kam- 
LAND is on the order of 100 events per year [29] , 
uniformly distributed across the detector. At SNO+, 
the reactor background would be a factor of about five 
lower [30 . There is 9.4 m of passive and active shield- 
ing in between the end of the sleeve and the beginning 
of the fiducial volume, including an instrumented water 
veto detector. This shielding is adequate for attenuat- 
ing/eliminating beam-related neutrons that can produce 
an IBD-like background, especially in consideration of 
the 3 MeV prompt energy threshold requirement. Fur- 
thermore, the sinusoidal- wave-like nature of an expected 
oscillation signal in L/E cannot be mimicked by back- 
ground. 

To perform an oscillation analysis, the antineutrino 
travel distance (L) and energy (E) must be recon- 
structed simultaneously on an event-by-event basis. 
Using KamLAND's detection capability as an exam- 



ple for the performance of a large scintillator detec- 
tor, the energy can be reconstructed with a reso- 
lution of 6.4%/s/E (MeV) [23]. With the antineu- 
trino event vertex in the detector known to within 
12 cm./ E (MeV) (23], the L resolution is dominated 
by the spatial distribution of activated 8 Li isotopes in- 
side the target and sleeve. The antineutrino creation 
point is distributed in the beam coordinate z according 
to an approximately uniform distribution, spanning the 
length of the 150 cm long sleeve. Although the spread in 
z dominates the smearing of the antineutrino baseline L, 
the distribution in terms of the transverse coordinates is 
taken into account as well. 

The IsoDAR oscillation analysis follows the method of 
Ref. [3T]. This analysis exploits the L/E dependence of 
oscillations, since L and E can be precisely reconstructed 
in the detector described. Eq. [T|is a good approximation 
for 3+1 (three active plus one sterile) disappearance fits 
to data, but in the case of 3+2 (three active plus two 
sterile), Eq. [T]is modified to accommodate Am 2 1; Amf^, 
and Am| 5 oscillations. If U is the mixing matrix, then 
the disappearance probabilities in the 3+1 and 3+2 sce- 
narios are given by 

P 3+ i = l-4|C/ e4 | 2 (l-|[/e4| 2 )sin 2 (Am 2 1 L/i;) (2) 
^3+2 = l-4[(l-|L/ e4 | 2 -|{/ e5 | 2 ) x 

(\U e4 \ 2 sm 2 (Aml 1 L/E) + 

\U e5 \ 2 sin 2 (Am 2 1 L/E)) + 

\U ei \ 2 \U e5 \ 2 sin 2 (Am 2 A L/E)]. (3) 

This assumes that contributions to disappearance from 
the [i and r elements of the mixing matrix U^, 
U T 4, and Urs) are negligible. We note that the current 
global fit improves significantly in the case of two sterile 
neutrinos [32 . 

Figure [3] illustrates the L/i?-dependent signal for ex- 
ample 3+1 and 3+2 oscillation signals after five years 
of running. The observation of an oscillation wave, fea- 
turing multiple peaks and valleys for currently favored 
values of Am 2 , makes this a highly compelling analy- 
sis. The wave also allows differentiation between 3+1 
and 3+2 models in most oscillation scenarios. The 3+2 
model-based oscillation probability shown in Fig. [3] uti- 
lizes the oscillation parameters given in Ref. [32]. These 
parameters represent the best fit of the world's appear- 
ance and disappearance data. 

IsoDAR can quickly probe the oscillation parameter 
space indicative of one or more sterile neutrinos. As the 
antineutrino source described can be constructed within 
five years, we compare the IsoDAR 95% CL sensitivity 
to experiments that can be accomplished on this time 
scale. The global fit region, encompassing all appear- 
ance and disappearance measurements, is shown along 
with this comparison in Fig. [4] Note that the global 
fit |32j pulls the reactor anomaly allowed region signifi- 
cantly lower in Am 2 due to the LSND and MiniBooNE 
appearance results, resulting in a Am 2 ~l-2 eV 2 . As can 
be seen in Fig. [4j the statistics-limited IsoDAR sensitiv- 



4 




FIG. 3: The L/E dependence of two example oscillation signatures after five years of IsoDAR running. The solid curve is 
the oscillation probability with no smearing in the reconstructed position and energy. The 3+2 example (right) represents 
oscillations with the global best fit 3+2 parameters from Ref. [32] , 



ity covers the 3+1 allowed range [sin 2 26 = 0.067 and 
Am 2 = 1 eV 2 ] at 20a in five years of running. Iso- 
DAR can rule out sin 2 20=0.067, Am 2 = 1 eV 2 at ha 
in 4 months. The "shape-only" and "rate+shape" based 
sensitivities are shown in the plot. It is clear that the flux 
normalization uncertainty is only relevant for oscillation 
sensitivity at high Am 2 (>15 eV 2 ), a region where the 
rapid oscillation wave becomes difficult to resolve. 
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trino oscillations near Am 2 ~ 1 eV 2 . The experiment can 
perform compelling tests of models for new physics that 
explain high Am 2 oscillations through the introduction 
of one or more sterile neutrinos. In a 3+1 model, Iso- 
DAR can address the global fit region for electron flavor 
disappearance at 20<r (5a) in five years (four months). In 
addition, the form of the oscillation wave can be recon- 
structed, allowing differentiation between the existence 
of 3+1 and 3+2 neutrinos. The large event sample also 
provides the possibility of a wide variety of other stan- 
dard model tests, including a precise measurement of the 
weak mixing angle. 



FIG. 4: The sensitivity of the IsoDAR experiment to elec- 
tron antineutrino disappearance in a five year physics run. 
The sensitivities for both rate+shape (solid line) and shape- 
only (dashed line) are shown. The ^iDAR [33] exclusion curve 
and Reactor+Gallium |3i allowed region are also shown, along 
with the expected sensitivities from the PBq source [34] and 
KATRIN [3H] experiments. 



The IsoDAR technique provides a high-intensity, low- 
energy source of antineutrinos with sensitivity to antineu- 
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